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Abstract 

In this Letter, we consider the Universe at the late stage of its evolution and 
deep inside the cell of uniformity. At these scales, the Universe is filled with in- 
homogcncously distributed discrete structures (galaxies, groups and clusters of 
galaxies). Supposing that the Universe contains also the cosmological constant 
and quintessence/phantom fluids with the constant equation of state parame- 
ter to, we investigate scalar perturbations of the FRW metrics due to inhomo- 
geneities. Our analysis shows that, to be compatible with the theory of scalar 
perturbations, the quintessence/phantom, first, should be clustered and, sec- 
ond, should have the equation of state parameter uj = —1/3. Therefore, this 
quintessence neither accelerates nor decelerates the Universe. We also obtain the 
equation for the nonrelativistic gravitational potential created by a system of in- 
homogeneities. Due to the quintessence, the physically reasonable solutions take 
place for flat, open and closed Universes. Quintessence is concentrated around 
the inhomogeneities and results in screening of the gravitational potential. 

Keywords: quintessence, phantom, cosmological constant, late time 
acceleration, inhomogeneous Universe, scalar perturbations 



1. Introduction 

The accelerated expansion of the Universe at late stages of its evolution, 
found little more than a decade ago [H, Q , is one of the most intriguing puzzles 
of modern physics and cosmology. Recognition of this fact was the awarding of 
the Nobel Prize in 2011 to Saul Perlmuttcr, Adam Riess and Brian Schmidt. 
After their discovery, there were numerous attempts to explain the nature of 
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such acceleration. Unfortunately, there is no satisfactory explanation so far 
(see, e.g., the state of art in Q). According to the recent observations 0, [Bj], 
the ACDM model is the preferable one. Here, the accelerated expansion is 
due to the cosmological constant. However, there is a number of problems 
associated with the cosmological constant. Maybe, one of the main of them 
consists in the adjustment mechanism which could compensate originally huge 
vacuum energy down to the cosmologically acceptable value and to solve the 
coincidence problem of close magnitudes of the non-compensated remnants of 
vacuum energy and the energy density of the Universe at the present time 0. 
On the other hand, it is well known that perfect fluids with the equation of 
state parameters w < —1/3 can also cause the accelerated expansion of the 
Universe. Such fluids are called quintessence 

0, H H and phantom [13, Ell for 
— 1 < lu < and oj < — 1, respectively. Therefore, they can be an alternative to 
the cosmological constant. It is of interest to investigate the viability of these 
models. 

In our Letter, we consider the compatibility of these models with scalar per- 
turbations of the Friedmann-Robertson- Walker (FRW) metrics at late stages of 
the Universe evolution. We explore a simplified model with a constant equation 
of state parameter u = const. We demonstrate that, first, these fluids must be 
clustered (i.e. inhomogeneous) and, second, u> = —1/3 is the only parameter 
which is compatible with the theory of scalar perturbations. Therefore, this 
quintessence neither accelerates nor decelerates the Universe. We also obtain 
formulas for the nonrelativistic gravitational potential created by a system of 
inhomogeneities (galaxies, groups and clusters of galaxies). We show that due 
to the quintessence, the physically reasonable expressions take place for flat, 
open and closed Universes. If quintessence is absent, the hyperbolic space is 
preferred fl2| . Hence, quintessence can play an important role. 

The Letter is structured as follows. In Sec. 2, we consider scalar pertur- 
bations in the Friedmann Universe filled with the cosmological constant, pres- 
sureless dustlike matter (baryon and dark matter) and quintessence/phantom 
fluids. Here, we get the equation for the nonrelativistic gravitational poten- 
tial. In Sec. 3, we find solutions of this equation for an arbitrary system of 
inhomogeneities for flat, open and closed Universes. These solutions have the 
Newtonian limit in the vicinity of inhomogeneities and are finite at any point 
outside inhomogeneities. The main results are summarized in concluding Sec. 
4. 



2. Scalar perturbations of FRW Universe 

Homogeneous background. 

To start with, we consider a homogeneous isotropic Universe described by 
the FRW metrics 

ds 2 = a 2 (dry 2 - ~f a pdx a dx f3 ) = a 2 (dry 2 - d X 2 - Ys 2 { X )d9 2 ) , (1) 
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where 

!sinx, xG[0,tt] for K = +1 

X, xe[0,+oo) for K = (2) 

sinhx, xG[0,+oo) for K, = — 1 

and /C = — 1,0, +1 for open, flat and closed Universes, respectively. As matter 
sources, we consider the cosmological constant A, pressureless dustlikc matter 
(in accordance with the current observations [H, [B| , we assume that dark matter 
(DM) gives the main contribution to this matter) and an additional perfect fluid 
with the equation of state p = cue where — 1 < uj < for the quintessence and 
uj < — 1 for the phantom fluid. In the present Letter, u> = const. The overline 
denotes homogeneous quintessence/phantom fluids. It can be easily seen from 
the conservation equation that in the case of the homogeneous perfect fluid 

3(l+w) 

I=£ \l(l+u>) > ( 3 ) 

where ao is the scale factor at the present time and eo is the current value of 
the energy density e. 

Because we consider the late stages of the Universe evolution, we neglect the 
contribution of radiation. Therefore, the Friedmann equations read 

— i 5 '- = kT q + A + ke (4) 

and 

2H' + H 2 +IC , 

5 = A — kuje , (5) 

a z 

where T~L = a' /a = (da/dr})/a and k = SttGn/c 4 " (c is the speed of light and 



Gat is the Newton's gravitational constant). Here, T is the energy-momentum 
tensor of the average pressureless dustlike matter. For such matter, the energy 
density T = ~pc 2 /a 3 is the only nonzero component, p = const is the average 
rest mass density [l2j]. It is worth noting that in the case IC = the comoving 
coordinates x a may have a dimension of length, but then the conformal factor 
a is dimcnsionlcss, and vice versa. However, in the cases /C = ±1 the dimension 
of a is fixed. Here, a has a dimension of length and x a are dimcnsionlcss. For 
consistency, we shall follow this definition for JC = as well. For such choice of 
the dimension of a, the rest mass density has a dimension of mass. 

It is well known that quintessence can provide the accelerated expansion of 
the Universe if the equation of state parameter satisfies the condition uj < —1/3. 
We can easily see it if we rewrite equations (0| and ([5]) as follows: 



-„4 A „2 „.„2 „ „3(l+w) 

(l + 3w), (6) 



a npc Ac kc £o<2 



a 6a 3 3 6 a 3 ^") 

where overdots denote the differentiation with respect to synchronous time t 
connected with conformal time rj: cdt = adr\. 
Scalar perturbations. 
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Obviously, the inhomogeneities in the Universe result in scalar perturbations 
of the metrics ([TJ. In the conformal Newtonian gauge, such perturbed metrics 



is 



IlEIl 



ds 2 w a 2 [(1 + 2<P)drf - (1 - 2^)j af )dx a dxP] , (7) 

where scalar perturbations $,\E f <C 1. Following the standard argumentation, 
we can put $ = vp. We consider the Universe at the late stage of its evolution 
when the peculiar velocities of inhomogeneities (both for dustlike matter and 
quintessence/phantom fluids) are much less than the speed of light: 

dx a dx a 1 v a 

— =a— - = — «1. 8 

drj dt c c 

We should stress that smallness of the nonrelativistic gravitational potential $ 
and peculiar velocities v a arc two independent conditions (e.g., for very light 
rclativistic masses the gravitational potential can still remain small). Under 
these conditions, the gravitational potential $ satisfies the following system of 
equations (see [12j for details): 

A$ - 3"H($' + U&) + 3/C$ = ^na 2 5T° + -m 2 5e , (9) 

JL(<S>> + H<S>)=0, (10) 

$" + 3U& + (2H' + H 2 )<S> - /C$ = ^na 2 5p , (11) 
where the Laplace operator 

*=^M^<"^) (i2) 

and 7 is the determinant of 7 Q| g . Following the reasoning of [l2j , we took into 
account that peculiar velocities of inhomogeneities are nonrelativistic and the 
contribution of 6Tp is negligible compared to that of ST® both for dustlike 
matter and quintessence/phantom. In other words, account of 8T® is beyond 
the accuracy of the model. This approach is completely consistent with fl5| 
where it is shown that the nonrelativistic gravitational potential is defined by 
the positions of the inhomogeneities but not by their velocities (see Eq. (106.11) 
in this book). In the case of an arbitrary number of dimensions, a similar 
result was obtained in [3- On the other hand, the motion of nonrelativistic 



inhomogeneities is defined by the gravitational potential (see, e.g., [17j). The 
perturbed DM remains nonrelativistic (prcssureless) that results in the condition 
STp = 0. For the quintessence/phantom we have 8T^ = —SpSp. In Eq. ©, 
Se is a fluctuation of the energy density for quintessence/phantom, while 5Tq is 
related to the fluctuation of the energy density of dustlike matter and has the 
form 0: 
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where Sp is the difference between real and average rest mass densities: 



8p = p-p. (14) 

From Eq. (fT0|) we get 

$(„,r) = -|£L, (15) 

where ip(r) is a function of all spatial coordinates and we have introduced c 2 
in the denominator for convenience. Below, we shall see that tp(v) ~ 1/r in 
the vicinity of an inhomogcncity, and the nonrclativistic gravitational potential 
$(77, r) ~ l/{ar) = 1/R, where R — ar is the physical distance. Hence, $ has the 
correct Newtonian limit near the inhomogeneities. Substituting the expression 
(fT5"|) into Eqs. © and (fTTj) . we get the following system of equations: 

i (A<p + 3K.tp) = ^nc 2 6T° + ^nc 2 Se , (16) 

\(H' - H 2 - K.)<p = \nc 2 5p. (17) 
a 6 2 

From the Friedmann equations (jj) and we obtain 

\ (W -n 2 -K) = — {-ijf a -k(1 + u)s) . (18) 
Then, Eq. (fTTJ) reads 

-k^ - «(1 + ^^j) V = nc 2 u5e . (19) 
a 4 ao a 4+JLJ / 

It should be noted that we consider quintessence/phantom fluids without ther- 
mal coupling to any other type of matter. It means, in particular, that evolution 
of its homogeneous background as well as scalar perturbations occurs adiabati- 
cally or, in other words, without change of entropy. Therefore, we preserve the 
same linear equation of state Sp = uSe with the same constant parameter uj for 
the scalar perturbations Sp and Se of pressure and energy density respectively, 
as for their background values p and e. 

Taking into account the expression (|15p . we get that in the right hand side 
of Eq. (p~3|) the second term is proportional to 1/a 4 and should be dropped 
because we consider the nonrclativistic matter. This is the accuracy of our 
approach, i.e. for the terms of the form of l/a ra , we drop ones with n > 4 
and leave terms with n < 4. Obviously, 4 + 3u < 4 for u < 0. Hence, we can 
draw the important conclusion regarding the purely homogeneous non-clustered 
quintessence/phantom fluids with 5p,Se = 0. For these fluids, we arrive at a 
contradiction because in Eq. (fl~9| the right hand side is equal to zero while the 
left hand side is nonzero. Therefore, such fluids are forbidden. Quintessence 
and phantom should be capable of clustering. In the papers (j| [l8| , it was also 
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pointed out that the quintessence has to be inhomogeneous. For inhomogeneous 
quintessence/phantom we get from Eq. (fl~9)) that 



5£ = —^^\^V- (20) 

It is worth noting that in our model neither the nonrelativistic gravitational 
potential $ ~ \ ja nor the quintessence/phantom density contrast Se/e ~ 1/a 
diverge with time (with the scale factor a) in spite of the negative sign of the 
ratio 5p/5e which is often treated as the speed of sound squared. Substituting 
(|2"0|) into (fH>)) . we obtain within our accuracy 

1 , A ± „ r s 1 2<^ 2 1 2 l+U 

r ( A(p + 3/Cw) = — kc — ^ —kc — ^ — £o , ; o y 

a 3 v ^ ^ ; 2 a 3 2 c 2 w a 4+3w ^ 

1 1 + u) a 1+3uj 
=> Aip + 31C(p = -nc A 5p — K£o al ° +3u> ip. (21) 

In this equation, all terms except the last one do not depend on time. Therefore, 
uj = — 1/3 is the only possibility to avoid this problem. Hence, we arrive at the 
following important conclusion. At the late stage of the Universe evolution, 
quintessence/phantom fluids are compatible with the scalar perturbations only 
if, first, they are inhomogeneous, and, second, they have the equation of state 
parameter cu = —1/3. Such fluid is quintessence which cannot provide the late 
time acceleration of the Universe (i.e. this quintessence neither accelerates nor 
decelerates the Universe (see Eq. ©)). 

For uj = —1/3, the equation for the gravitational potential and the fluctua- 
tion of the energy density of the quintessence read, respectively: 

+ (siC - ^^eoalj <p - 4ttG w (p-p) (22) 

and 

In the next section, we shall investigate Eq. (j2"2"j) depending on the curvature 
parameter JC. We shall show that reasonable expressions of the conformal gravi- 
tational potential <p exist for any sign of JC. This takes place due to the presence 
of the quintessence with u = —1/3. If quintessence is absent, the hyperbolic 
model fC = — 1 is preferred [l2|. Therefore, the positive role of quintessence is 
that its presence gives a possibility to consider models for any JC. 

3. Gravitational potentials 

It is convenient to rewrite Eq. (f22|) as follows: 

A(j)- X 2 4> = AttGnp, (24) 
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where the truncated gravitational potential is 



* = *-^, A^O, (25) 



and 

2 _ 8ttG n 2 
A = — j — e a - 3/C . (26) 



On scales smaller than the cell of uniformity size (which is of the order of 150 
Mpc) and on late stages of evolution, the Universe is filled with inhomogeneously 
distributed discrete structures (galaxies, groups and clusters of galaxies) with 
dark matter concentrated around these structures. Then, the rest mass density 



p reads [12 1 



where m, is the mass of i— th inhomogeneity. Therefore, Eq. ([Ml) satisfies the 
very important principle of superposition. It is sufficient to solve this equation 
for one gravitating mass nii and obtain its gravitational potential fa. The 
gravitational potential for all system of inhomogencities is equal to a sum of 
potentials fa. It is worth recalling that the operator A is defined by Eq. (|T2j) . 
As boundary conditions, we demand that, first, the gravitational potential of 
a gravitating mass should have the Newtonian limit near this inhomogeneity 
fa ~ \ jr and, second, this potential should converge at any point of the Universe 
(except the gravitating mass position). 

It seems reasonable to assume also that the total gravitational potential 
averaged over the whole Universe is equal to zero (see, e.g., [Hj]): 



- AttG n p — yr^ 1 



A 2 ' Y 



, <P = J2 v / & dV > ( 28 ) 



where V is the volume of the Universe. This demand results in another physi- 
cally reasonable condition: 8e = (see Eq. ([23]) "). 



Flat space: JC — 0. 



In the case Eq > — > A 2 = S ^ N £o a o > 0, the solution of (|24|) for a separate 
mass 171% satisfying the mentioned above boundary conditions reads 

= nttm exp ^_ Ar ) A>0,0<r<+oo. (29) 

r 

It can be easily seen that this truncated potential has the Newtonian limit 
for r — > 0. This expression shows that quintessence results in the screening 
of the Newtonian potential. A similar effect for the Coulomb potential takes 
place in plasma. In our case, the screening originates due to specific nature 
of quintessence. It is worth mentioning that the exponential screening of the 
gravitational potential was introduced " by hand" in a number of models to solve 
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the famous Seeliger paradox (see, e.g., the review [19|). In our model, we resolve 
this paradox in a natural way due to the presence of quintessence. 

For a many-particle system, the total gravitational potential takes the form 



-G N V 8 exp (-A|r - r t \) 
* — ' r — r,i 



A 2 



(30) 



Substituting ([3"0)l into ([25]) , we get for the fluctuation of the quintessence energy 
density the following expression: 



Se = 



■exp(-A|r-rj|) 



c 4 P 
2e Q a 2 l 



(31) 



Therefore, we arrive at a physically reasonable conclusion that these fluctua- 
tions are concentrated around the matter/dark matter inhomogeneities and the 
corresponding profile is given by Eq. ()31l) . 

The averaged value of the i-th component of the truncated potential over 
some finite volume Vq is 



4-7T 

Vq 



ii 



Gnth. 



exp(— Xr) 



r 2 dr 



(32) 



V, 



(i 



cxp(-Arp) 
A 



ro 



Then, letting the volume go to infinity (ro 
gravitating masses, we obtain 

<p = -Gn~p 



1 

" A 

Vq 



1 

A 2 



4tt 

A 2 ' 



-oo ) and taking all 
(33) 



where p = lim ^rrii/Vo. Therefore, the averaged gravitational potential 

Vb-f+oo a 



([28]) is equal to zero: Tp = 0. Consequently, 5e = 0. 

The case £o < =>■ A 2 = — /j 2 < is not of interest. Here, we get the 
expression (f>i = — {GNTn.i/r) cos(/ir) which does not have clear physical sense. 
Additionally, this expression does not allow the procedure of averaging. 



Spherical space: K, 



-1. 



Let us consider, first, the case A 2 = Sn ^l N eoa 2 , — 3 = — p? < 0. This case is 
of interest because it allows us to perform the transition to small values of the 
energy density of quintessence: eo — > 0. Here, the solution of ([24]) for a separate 
mass m,; is 



= —Gjsrrrii- 



sin 




+ 1 


sin 




sinx 



< % < 7T. 



(34) 
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For i//i 2 + 1 7^ 2, 3, . . . (we would remind that /i 2 ^ 0), this formula is finite 
at any point % G (0, tt] and has the Newtonian limit for x — > 0. In the case of 
absence of quintessence e = — > ^//i 2 + 1 = 2, this expression is divergent at 



X = 7r. We demonstrated this fact in our paper [12(. Therefore, quintessence 
gives a possibility to avoid this problem for the models with JC = +1. It can 
be easily verified that for the total system of gravitating masses, the averaged 
value of the total truncated potential has the form of (|3"3")l that results in Tp = 

In the case A 2 > 0, the formulas can be easily found from (|34|) with the help 
of analytical continuation fx —> In other words, it is sufficient in Eq. (|34[) to 
replace /j, 2 by —A 2 . The obtained expression is finite for all x € (0, 7r] and the 
averaged gravitational potential is equal to zero: Tp = Se = 0. 

Hyperbolic space: K, = — 1. 

Here, the most interesting case corresponds to A 2 = eoa 2 + 3 > 0. This 
choice of sign gives a possibility to perform the transition to small values of the 
energy density of quintessence: £o — > 0. Then, the desired solution of Eq. (|24|) 
for a mass m, is 



- . r exp 
sinn x 



(-XV 7 A 2 + l) , 0<x<+^. (35) 



If q uintessence is absent (eo = 0), then we reproduce the formula obtained in 
[12j. On the other hand, the expression (f3"5l) shows that for £o >0^A 2 + 1 > 4, 
quintessence enhances the screening of the gravitating mass. For a many-particle 
system, the total gravitational potential takes the form 



cxp(-Z 4 v / A r +I) , 4ttG n p 
l p = -G N ^m i — + ^^, (36) 

i 

where li denotes the geodesic distance between the z-th mass rrii and the point 
of observation. Similarly, using Eq. Q34p. we can write the expression for the 
total potential in the case of the spherical space. 

Taking into account that the averaged total truncated potential has again 
the form (|3"3"|) . the procedure of averaging leads to the physically reasonable 
result: Tp — =>■ Se = 0. 

Concerning the case A 2 < 0, the truncated gravitational potential is finite 
in the limit x —> +oo. However, the procedure of averaging does not exist here. 
Therefore, this case is not of interest for us. 

To conclude this section, we discuss briefly the case A 2 = 0. For K = 0, —1, 
the principle of superposition is absent now. To make the gravitational potential 
finite at any point including the spatial infinity, we need to cutoff it smoothly at 
some distances from each gravitating mass. If K, = 0, then quintessence is absent 



and this case was described in detail in [12(. It was shown that the averaged 



gravitational potential is not equal to zero. This is a disadvantage of such 
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models. In the case K, = +1, the principle of superposition can be introduced 
due to the finiteness of the total volume of the Universe. Here, the comoving 
averaged rest mass density can be split as follows: p = TOj/(27t 2 ) = J2i Pi- 
Then, Eq. (|22j) can be solved separately for each combination (mi,/5j). As a 
result, the gravitational potential of the z-th mass is 

G N mi cos X / -, X\ n . s ( o-7\ 

tpi = — G N rrii- 1 , 0<x<tt- (37) 

2n sinx ^ 7r/ 

This potential is convergent at any point x 7^ 0, including x = ^- It is not 
difficult to see that Tp i = 0. Therefore, the total averaged gravitational potential 
is also equal to zero: Tp = £\ Tp i = & = 0. 

4. Conclusion 

In our Letter, we have investigated the role of quintessence and phantom 
field for the Universe at late stages of its evolution. It is well known that these 
fields can be an alternative to the cosmological constant explaining the late time 
acceleration of the Universe. It happens if their parameter of the equation of 
state u> < —1/3. To check the compatibility of these fluids with observations, 
we consider our Universe at scales much less than the cell of homogeneity size 
which is approximately 150 Mpc. At such distances, our Universe is highly 
inhomogeneous and the averaged Friedmann approach does not work here. We 
need to take into account the inhomogencities in the form of galaxies, groups 
and clusters of galaxies. All of them perturb the FRW metrics. 

To clarify the role of quintessence and phantom, we endowed the Universe 
with these fluids with a constant equation of state parameter u. We have 
shown that quintessence and phantom are compatible with the theory of scalar 
perturbations if they satisfy two conditions. First, these fluids must be clustered 
(i.e. inhomogeneous). Second, the parameter of the equation of state w should 
be —1/3. Therefore, this quintessence neither accelerates nor decelerates the 
Universe. 

Then, we have obtained the equation for the nonrelativistic gravitational po- 
tential. We have shown that due to quintessence the physically reasonable solu- 
tions take place for flat, open and closed Universes. The presence of quintessence 
helps to resolve the Seeliger paradox [l9[ for any sign of the spatial curvature 
parameter K,. If quintessence is absent, the hyperbolic space is preferred [l2j |. 
Hence, quintessence can play an important role. Quintessence is concentrated 
around the inhomogeneities and results in screening of the gravitational poten- 
tial. 
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